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%-( ■ Abstract 
Oh' 

' The paper presents an ant colony optimization metaheuristic for collaborative planning. Collabora- 

tive planning is used to coordinate individual plans of self-interested decision makers with private 
information in order to increase the overall benefit of the coalition. The method consists of a 
new search graph based on encoded solutions. Distributed and private information is integrated 
via voting mechanisms and via a simple but effective collaborative local search procedure. The 
approach is applied to a distributed variant of the multi-level lot-sizing problem and evaluated by 
C/3 ' means of 352 benchmark instances from the literature. The proposed approach clearly outperforms 

^ existing approaches on the sets of medium and large sized instances. While the best method in 

the literature so far achieves an average deviation from the best known non-distributed solutions 
of 46 percent for the set of the largest instances, for example, the presented approach reduces the 
' average deviation to only 5 percent. 

j. ■ Keywords: Group decisions and negotiations; metaheuristics; collaborative planning; ant colony 

, optimization; lot sizing 
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1 Introduction and literature review 

A coalition of decision makers with private information that requires collaborative production planning 
is considered. Each of the decision makers, hereafter referred to as agents, is selfish and seeks to 
. implement his or her local optimal production plan to minimize his or her local costs. However, the 
agents could improve their situation by coordinating their individual production plans. The cost savings 
due to a better global production plan could be allocated to the agents in order to overcompensate 
them for deviating from their local optimal production plans. From a more general point of view, the 
problem of finding an appropriate allocation of cooperative gains is dealt with by cooperative game 
theory. The special situation of joint planning of selfish agents in case of private information considered 
here is the subject of collaborative operations planning (Dudek and Stadtler, 2005). 

The state of the art in collaborative planning is discussed by Frayret (2009) and by Stadtler 
(2009). Overviews of automated negotiation approaches relevant to collaborative planning are given 
by Lomuscio et al. (2003) and Strobel and Weinhardt (2003). Applications of collaborative planning 
to supply chains are discussed by Ertogral and Wu (2000), and Fink (2004, 2006), for example. 

The proposed collaborative planning approach is applied to a distributed variant of the well-known 
multi level uncapacitated lot sizing problem (MLULSP) introduced by Yelle (1979). The distributed 
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MLULSP (DMLUSLSP) was presented by Homberger (2010) and assumes several local and selfish 
agents with private information instead of a single agent with full information. The DMLULSP cov- 
ers some important features of real world problems. There are several final products, a multi level 
production structure, and a trade-off between inventory and setup costs. Coordination is difficult, 
because there are agents with private information and conflicting objectives. Finally, the problem is 
computationally challenging, as it is NP-hard for general product structures (Arkin et al., 1989). As 
to computational experiments, a variety of benchmark instances is available in the literature which 
allow for comparative test including other approaches. 

Lot Sizing Problems are of high relevance in modern supply chains (VoEand Woodruff, 2006). 
Surveys of the rich literature on lot sizing problems are given for example by Drexl and Kimms 
(1997), Jans and Degraeve (2007), or Winands et al. (2011). Integration of lot sizing with cooperative 
game theory is discussed by Drechsel and Kimms (2010, 2011). Metaheuristic solution approaches 
for the non-distributed MLULSP are amongst others presented by Dellaert and Jeunet (2000, 2003), 
Pitakaso et al. (2007), Homberger (2008), Xiao et al. (2011), and Xiao et al. (2012). 

The literature discusses at least two distributed lot sizing variants. The one studied by Lee and Kumara 
(2007b, a) is based on auctions and has been tested on small size instances. For the other, the distributed 
MLULSP, several metaheuristic approaches that integrate concepts from group decision making are 
known in the literature: a simulated annealing approach (SA, Homberger, 2010), an approach based 
on an evolutionary strategy (ES, Homberger, 2011) and one based on ant colony optimization (AGO, 
Homberger and Gehring, 2010, 2011). These approaches differ by the used metaheuristic principle, 
the applied negotiation or voting mechanisms, and the relevant objective functions (minimizing total 
global cost or maximizing fairness). ES is especially useful for maximizing fairness. SA and AGO 
strive to minimize total global cost. In terms of solution quality, AGO outperforms SA on medium 
sized instances (40 and 50 items with 12 and 24 periods) but not on small or large instances (500 items, 
and 36 or 52 periods). Furthermore, AGO is significantly slower than SA on small, medium, and large 
instances. 

The aim of this paper is to present a new collaborative ant colony metaheuristic (GAGM) for solving 
the DMLULSP under global total cost minimization. Compared to Homberger and Gehring (2011) and 
the other mentioned approaches, the proposed GAGM incorporates the following new conceptual and 
methodical ideas: 

1. Encoded solutions are represented by a new and simplified search graph. Let m denote the 
number of items to produce and let n denote the number of periods. The search graph proposed 
by Homberger and Gehring (2009) and used in Homberger and Gehring (2011) requires m(n + l) 
nodes and ^mn(n + 1) directed edges while the new search graph uses 2mn + 1 nodes but only 
4mn — 2 edges. The new search graph facilitates learning of pheromone values and therefore 
makes it easier to discover promising paths which lead to high quality solutions. 

2. Individual preferences of the agents - which are private - are aggregated by a voting mechanism 
and by a new collaborative local search procedure. Usually, to be able to reach a (near) optimal 
solution fast, an ant colony metaheuristic requires heuristic information on the problem, for 
example, the costs of the agents. However, this operational information is sensitive and therefore 
private. The new approach exploits this information indirectly, which enables to guide the search 
via heuristic information and, nevertheless, keeps operational information private. 

3. The new approach does without rank based voting mechanisms. This reduces computational ef- 
fort and increases robustness, that is, there are less opportunities for undesired strategic behavior 
of the agents. 

As it is shown below, the incorporation of these concepts has a high positive impact on the solution 
quality of the resulting approach. 

The rest of the paper is organized as follows. In Section 2 the considered distributed multi- level lot 
sizing problem (DMLULSP) is characterized. Section 3 describes the developed collaborative ant colony 
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Table 1: Notation for the DMLULSP. 



Problem parameters 



M 


p\ mimfiPTitlv IjirP'P TiiiTnhpT 


Ul 


iiuniber of items 


n 


number of production periods 


I 


set of items, / = {1, . . . , m} 


T 


set of possible production periods, T = {1, . . . , n} 


Si 


setup costs per period for item z S / 


hi 


inventory holding costs per period and per unit of item z G / 


U 


lead time required to assemble, manufacture, or purchase item i S J 




number of items i required to produce one unit of item j with i, j E I , i ^ j 


r+{i) c / 


set of all direct successors of item z G / 


r-(z)c/ 


set of all direct predecessors of item i E I 




exogenous demand (unit of quantity) for item i E {j E /|r+(j) = 0} in 




period t eT 


Decision variables 




endogenous demand (unit of quantity) of item i E {j E I\T~^{j) ^ 0} in 




period t eT 


lit 


inventory (unit of quantity) of item z G / at the end of period t E T 


Xit 


lot size (unit of quantity) of item z G / in period t E T 


Vit 


binary setup decision, yit = 1 if item z G / is produced in period t E T and 




Uit = otherwise 



metaheuristic for solving the DMLULSP. Subject to Section 3 is the evaluation of the new approach 
by means of a computational study with 352 benchmark instances. Finally, concluding remarks are 
given in Section 5. 

2 A distributed multi-level uncapacitated lot-sizing problem 
2.1 Classical centralized problem formulation 

The distributed lot-sizing model studied in this paper extends the well-known multi-level uncapacitated 
lot sizing problem (MLULSP, cf. Yelle 1979, Steinberg and Napier 1980, Dellaert and Jeunet 2000). 
Therefore, the MLULSP is presented first using the notation given in Table 1. The MLULSP assumes 
a single decision maker who is aware of all information required for planning, especially, the setup costs 
and inventory holding costs per item and period. A formal description of the MLULSP is given by the 
formulas (1) to (8). 



min /"'^(y) = XI " + ^» ' 




(1) 








S. t. lit = ^i,t-l + Xit — dit , 


vi £i,yt£ T, 


(2) 


k,o = , 


Vi G /, 


(3) 


lit > , 


ViG/, vtGr\{o}, 


(4) 


dit = X rij ■ Xj^t+ti , 


Vi G {i G / 1 r+(j) / 0}, t G T, 


(5) 








Xit -M-yit<0, 


yi £i,yte T, 


(6) 


Xit>0, 


Vi £l,yte T, 


(7) 


Vit G {0, 1} , 


Vz G /, Vt G T. 


(8) 
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Figure 1: Example of a product structure and a partition of the items to two agents (cf. Homberger, 
2010). 

The goal of the MLULSP is to minimize the total costs /'^'^ of a single, central decision maker. 
These are expressed by the objective function (1) which sums up the setup costs and the stockholding 
costs for all items i G / over all periods t € T. The inventory balance is guaranteed by (2). For all 
items, the inventory of the first period t = is zero (3) and for remaining periods non- negative (4). 
For each period, the demands da for the level zero items z € / with r~''(i) = are given. The demands 
for the remaining items are determined by (5). These constraints ensure that the production of item j 
in period t + tj triggers a corresponding demand da for all i G r~(j), that is a demand for each item i 
that is preceding item j in the bill of materials. Without loss of generality, rjj- = 1 is assumed in (5). 
The lot size xa is non-negative (7). If xa > 0, that is, item i is produced in period t, then yit = 1, 
otherwise ya = 0. This is enforced by the constraints (6) and (8). 

Arkin et al. (1989) have shown that the MLULSP is NP-hard for general multi-level product struc- 
tures, i.e., for product structures, where each item can have multiple successors and predecessors like 
in Figure 1 (Salomon and Kuik, 1993). Figure 1 shows an example for a graphical representation of 
a product structure with m = 40 which follows the literature (Bookbinder and Koch, 1990). The 
production dependencies between items are indicated by arrows. Furthermore, items are classified by 
decreasing production level from top to bottom. Each end product f G / with T~^{i) = is located at 
the highest level which is level 0. Each common part is located at the lowest level at which it is used 
anywhere in the product structure (Dellaert and Jeunet, 2000). The allocation of items to a group of 
agents also shown in Figure 1 is described next. 

2.2 Extension to a distributed group decision problem formulation 

Following Homberger (2010), the MLULSP is now extended to a distributed group decision problem 
which is called distributed multi- level uncapacitated lot-sizing problem (DMLUSLSP). Instead of a 
single decision maker or agent who is responsible to produce all items /, the responsibility is jointly 
assigned to a group of agents A. Each agent a G A is responsible to produce the set of items with 
UaGA = I and PlaGA ~ ^'^^ instance, an agent might represent a profit center or an independent 
company and the agents might interact in a supply chain. 

All agents are autonomous and self-interested. Therefore, the individual objective function fa of 
agent a G A is to minimize his or her total local costs for producing the items /„, that is. 
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fa{y) = ^ "^(Si ■ Vit + hi ■ lit). 



(9) 



Furthermore, for all agents a ^ A and for all items i € la, the DMLULSP assumes asymmetric 
information regarding the cost parameters Si and hi. These cost parameters are private information. 
That is, at the initiation of planning they are only known to agent a and during planning agent a 
does not want to reveal these to other agents. If this information would be common knowledge, price 
negotiations between agents might be negatively affected. On the other hand, symmetric information 
is assumed regarding the bill of materials. This assumption can be justified by some kind of common 
industry knowledge or joint development of products (e.g. collaborative engineering). 

The DMLUSLSP consists of the constraints (2) to (8) and the objective function (10) which mini- 
mizes the total global costs 



In comparable scenarios, the total global costs are also considered of interest by Fink (2004), Zimmer 
(2004), Jung et al. (2005), and Dudek and Stadtler (2005, 2007). In the following, the function (10) 
is simply referred to as total cost function or central cost function, the function (9) is also denoted 
as individual or local cost function of agent a £ A. Usually, the minimization of the individual cost 
functions conflicts with the objective of minimizing the total costs. In order to resolve these conflicts 
and support the agents in agreeing on a joint production plan in order to produce all items /, the 
solution approach proposed in the following section uses collaborative local search and different voting 
mechanisms. 

3 A collaborative planning metaheuristic based on ant colony opti- 
mization 

3.1 Overview 

A metaheuristic solution approach based on ant colony optimization (AGO) is presented that espe- 
cially takes into account the distributed nature of the DMLULSP. The distributed solution approach 
is denoted as collaborative ant colony metaheuristic (CACM); an overview is given by Alg. 1. The 
metaheuristic CACM works on an encoded representation of a DMLULSP solution which is described 
in Sect. 3.2. Components of an encoded solution are represented by a search graph G which is traversed 
by artiflcial ants in order to construct encoded solutions. Sect. 3.3 introduces the search graph G and 
Sect. 3.4 presents the stochastic ant construction procedure. Each constructed solution is improved by 
a collaborative local search heuristic (cf. Sect. 3.5). Both, the used graph structure and the local search 
heuristic incorporate new concepts compared to the ant metaheuristic of Homberger and Gehring 
(2011). After a set B of encoded solutions has been constructed, the agents negotiate about the 
acceptance of a solution in B by casting votes according to one of several voting rules described in 
Sect. 3.6. Finally, the pheromone information related to the arcs of the search graph are updated 
according to the present accepted solution e* (Sect. 3.7). That is, in the next iteration of CACM ants 
use a different probability distribution in order to construct a new encoded solution. 

CACM requires several input parameters. The problem data include the number of items m € N, 
the number of periods n € N, and the set A of agents. The termination of CACM is controlled by the 
number of generated solutions s G N. The ballot size 6 € N, 6 4^ \B\, specifies the number of solutions 
which can be negotiated by the agents A in a single election. The remaining parameters influence the 
update of the pheromone information on the arcs of the search graph G. The minimum and maximum 
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Algorithm 1: Collaborative Ant Colony Metalieuristic 



end 

e* ^ voting(i?, A) 

T updatePheromoneValues(e*, r, t^' 
end 

return e* 



// Section 3.3 



Input: no. of items m, no. of periods n, set A of agents, no. of generated solutions s, ballot 
size b, min and max pheromone values r"**" and r™'^^, evaporation rate p, 
intensification rate a 

Output: jointly accepted encoded solution e* 

initialize search graph G and pheromone values r 
e*j ^ 0, i = 1, . . . ,171] t = 1, . . . ,n 
for s ^ 1 to s do 
B^{} 

for 6 1 to 6 do 

e ^ constructEncodedAntSolution(G, r) 
e ^ collaborativeLocalSearch(e, A) 

B ^ BLl{e} 
s s + 1 + m(n — 1) 



// Section 3.4 
// distributed, Section 3.5 



// distributed, Section 3.6 
%p,a) // Section 3.7 



pheromone values are defined by -j-™-™ ^ ^ r^j-^^ pheromone evaporation rate and the pheromone 
intensification are given by p and a with < p, o" < 1. 

The solution approach CACM can deal with asymmetric information distributed among the agents. 
CACM is executed by a neutral mediator. In order to construct encoded solutions, the mediator has to 
be aware of the items / to produce and of the possible production periods T. The agents do not have 
to reveal private information like costs or free production capacities, neither to the mediator nor to 
other agents. To control the search process in order to find a jointly accepted solution e* , the mediator 
has to interact with the agents in both steps marked as distributed in Alg. 1. CACM outputs a jointly 
accepted solution e*. 

3.2 Representation of a solution 
3.2.1 Solution encoding 

For the MLULSP, the decision variables xa represent the lot size of item i produced in period t. The 
binary decision variables yit represent the setup decisions which indicate wether item i is produced in 
period t at all. Because of the structure of the MLULSP, it is possible to determine the optimal lot 
sizes Xit, if the optimal setup decisions yu are known (Dellaert and Jeunet, 2000). 

Therefore, the proposed CACM focuses its search effort to approximate the set of optimal setup 
decisions. However, CACM does not operate on the y^ directly, but indirectly on an encoded repre- 
sentation of the setup decisions. For each binary setup decision variable yn, there is an encoded binary 
setup decision variable en following the suggestion of Homberger (2008). Both variables are related as 
follows: 



{en = 0)^ {yu = 0), (11) 

{eit = l)^{yit = 0\/yit = l). (12) 

Accordingly, the variables eu = 1 indicate a possible, but not mandatory, production of item z G / 
in period t & T. Another ACO approach that also uses encoded decision variables is suggested by 
Fink (2007), for example. Whether a possible production takes place if eu = 1, can be deduced by the 
following decoding rule. 
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Figure 2: Example of a search graph for two items and three periods. 
3.2.2 Solution decoding 

An encoded solution e € {0, l}™^" is transformed to a DMLULSP solution in a two steps procedure. 
The product level of an item i determines the order in which demand du , lot size xu , and inventory 
are calculated. The procedures starts with the lowest level items (final products) and advances with 
items of increasing product level. That is, the demands for the final products (level 0) are determined 
first, the demands for level 1 items are calculated next, and so on (cf. Homberger (2008)) 

Hence, the items are considered in the sequence of their non- decreasing product level. If item z is a 
final product (level 1), then da is given for all 1 < j < T. For each non- final product i the demand can 
be calculated according to equation (5) using to the demands of all successors of i indicated by T{i). 
All dkj,k G r(i) have already been calculated, because of the used calculation sequence. By means 
of the second decoding step, the setup decisions, lot sizes, and stored quantities are calculated. One 
should note, that with the used encoding and decoding rules, each encoded solution with a e^i = 1 for 
i = 1, . . . , N represents a feasible decoded solution in any case. 

3.3 Definition of a search graph 
3.3.1 Search graph 

The AGO principle is guided by the image of an ant that constructs a feasible solution by traversing a 
graph which consists of solution components specific to the problem. A solution is represented by the 
path chosen by the ant. hi order to enable such an artificial ant to construct an encoded solution for 
the DMLULSP, a directed graph G = {V, E) with a set V of nodes and a set E of directed edges is used 
(see Figure 2). The node set V consists of three disjunct subsets, i.e., V = {v} l^V^\J V^. The starting 
point of an ant is the initial node v. For each of the two possible values of an encoded decision variable 
eij there exists a node. The encoded decision to produce item i in period j {eij = 1) is represented by 
node vjj G . The encoded decision Sij = is represented by node vfj G V'^. Therefore, ClV^ = ^ 
and \V^\ = \V^\. The nodes in are called black nodes, the nodes in are called white nodes. 

The set of directed edges E is defined as follows. The initial node v is associated with two directed 
edges {v,vl i) and (v,v^-^). Furthermore, each node v ^ V \ {Vm m'^m n} ^^^s as the origin of exactly 
two directed edges. From a black node two edges originate, one points to another black node and one 
points to a white node; analogously, from a white node one edge points to a black node and one edge 
points to another white node. More precisely, let a,b € {0,1}, a ^ b and v"" G y",t>^ G V'^. The pair 
of edges originating from vfj is defined as {vfj,vfj_^_i) and (^i^-, t'^'j+i), provided that j < n. Else, if 
j = n, then the pair of edges is defined as {v^j,vfj^i i) and {vfj,v\j^i i). 

To ease presentation, two subsets of directed edges E^^ and E^ are introduced. E^ contains those 
edges that point to a black node and E^ contains those edges that point to a white node. 

E = E^UE^ with E^r\E^ = 0, (13) 
E^ ■.= {{i,j) (^E\i(^V ^3 (^V^}, (14) 
:= {(i,j) G G 1/ A j G (15) 
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An encoded solution to the DMLULSP is represented by all nodes v £ V UV on any path through 
G that starts at node v and ends either at node t^^.n '^^ node v^.n- 

3.3.2 Initialization of pheromone trail 

The sequence of directed edges traversed by an ant depends on the pheromone value T(^ij) placed on 
each edge G E. Higher pheromone values T(^ij) on an edge increase the probability that 

an ant traverses edge Homberger and Gehring (2009) proposed an AGO approach for a (non- 

distributed) MLULSP on a different search graph and successfully pursued the idea to initially construct 
each possible solution with equal probability ('balanced path'). In contrast to that idea, the goal of 
this approach is to initially construct unbalanced paths that predominantly consist of black nodes. 
Following the max-min ant system proposed by (Stiitzle and Hoos, 2000), a minimum and a maximum 
pheromone value for each edge (i, j) G E are defined as r™" and -7-™'"^. The initial pheromone value 
r(jj) for each edge € E^ is set to 1, and the initial pheromone value for G E^ is set to 

^max 

3.4 Construction of a solution 

3.4.1 Component selection probability 

In order to construct an encoded solution e, the graph G is traversed by an ant along the directed 
edges. The traversal of G starts at the initial node v. Until the ant arrives at node f^n or v^^, the 
ant decides randomly at each node v £ V \ {vj^ rp , v'j^ rp} which of the two possible edges {v,i) G E^ 
and {v,j) G E^ to choose. The random decision depends on the pheromone intensity T(j of an edge 
(i,j) G E. The probability to choose edge (i,j) is 

Pi^,J) = ^ ^ith T(,,,) + r(,,,) = r--^- and G E' A (i, k) £ E^ . (16) 

Hence, an according production for the first period of demand is planned and the calculation of a 
feasible solution is guaranteed (Dellaert and Jennet, 2000). 

3.4.2 Heuristic information 

Typically, AGO approaches from the literature additionally bias the selection probability of a compo- 
nent (here, an edge («, j) G E) via heuristic information. In the present case, however, the DMLULSP 
is a distributed problem with asymmetric information which is solved by a mediator. Because the 
mediator is neutral, he or she cannot take into account heuristic information that would depend on an 
individual agent and therefore privilege (or penalize) this agent. For this reason, equation (16) does 
without local heuristic information. 

3.5 Collaborative local search heuristic 

An encoded solution e constructed by an ant is improved by a collaborative local search heuristic (CLS, 
cf. Alg. 2). In this heuristic, a move is only executed, if each agent a £ A agrees. A move is defined 
as a bit shift, that is, for all i G / and all t G T, t > 1, the value of en is shifted from zero to one or 
one to zero, respectively. Each shift leads to new solution e' . The solution e' is submitted to all agents 
a £ A. If and only if e' has lower local costs fa{^') compared to /a(e) for each agent a G ^, the next 
move is executed on e' and so on. 

In CACM (cf. Alg. 1), the counter of generated solutions s is increased by l + m(n— 1), as there are 
m{n — 1) possible moves and each move applied to a generated ant solutions leads to a new solution. 
Furthermore, the final solution e generated by GLS is in GAGM again submitted to all agents. This 
time, the agents compare e to the jointly accepted solution e* and decide via voting, if e should become 
the new jointly accepted solution e*. 
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Algorithm 2: Collaborative local search 



Input: encoded solution e, set of agents A 
Output: possibly modified encoded solution e 

foreach t ^ 2, . . . ,n do 



foreach i 1, . . . , m do 




// bit shift 



s ^ s + l 
veto false 
foreach a ^ A do 



// complete approval voting 
// each agent a can veto e 



if /a(e') > fa{e) then veto true 



end 

if veto — false then e e' 
end 
end 

return e 



3.6 Distributed decision making by means of voting 

3.6.1 Complete approval voting 

In approval voting each agent can vote for as many candidate solutions in B as he or she wants 
and the solution with the most approvals wins (Brams and Fishburn, 1978). For the problem at 
hand, a modified and restricted version of approval voting is used which is denoted by the term 
complete approval voting. An encoded candidate solution e is only accepted as winner, if all agents 
a & A approve e. Therefore, by accepting e, no agent increases his or her current local costs C*, i.e., 
Ca(e) < C*, \/a £ A. For complete approval voting, the ballot size is restricted to 6 = 1. Although 
this voting rule is quite simple, it has not been used in the previous approaches of Homberger (2010, 
2011) and Homberger and Gehring (2011). 

3.6.2 Adapted Borda maximin voting rule 

Each agent a A ranks the candidate solutions e £ B according to his or her preferences. A ranking 
of e G -B by agent a is given by >- ep_i y . . . y ef where candidate solution is preferred most 
by agent a and receives p points. Solution ep_i is ranked second and receives p — 1 points, candidate 
solution is preferred least and receives 1 point. Given all rankings by all agents a £ A, the Borda 
maximin rule selects that candidate solution e' which maximizes the lowest number of points assigned 
by an agent to a candidate solution. 

Finally, the solution e' identified by Borda voting is compared to the jointly accepted solution e* by 
means of complete approval voting. The winning solution becomes the new jointly accepted solution 



3.6.3 Ravels or minimax voting rule 

The minimax voting rule, introduced by Rawls (1971), minimizes the maximum local cost of an agent. 
In contrast to the voting rules Borda maximin and complete approval, the Rawls voting rule requires 
the agents to reveal their local cost to the mediator. For each encoded candidate solution e £ B each 
agent a G A submits his or her local costs Ca{e) to the mediator. 



e*. 



e' = minmax{Ca(e)}. 



(17) 
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Finally, the solution e' identified by Rawls voting is compared to the jointly accepted solution e* by 
means of complete approval voting. The winning solution becomes the new jointly accepted solution 
e*. 

3.7 Pheromone update 

As the search advances, the pheromone information on each arc is continuously updated. First, the 
pheromone T(j j) of an edge S E evaporates by the evaporation rate p. Second, T^ij) is increased 
again by ar"^'^^, if and only if edge is part of the path chosen by an ant (that implies e^^^-^ = 1). 

Furthermore, equation (18) takes into account the smallest and largest possible pheromone values r'"*"' 
and r™""^, respectively. 

^ rmin{(l - p) ■ T(,j) + a ■ r'"'^-; r™-^-} if e^. = 1, 
^^'''^ ^ \max{(l - p) ■ mj); r™"} else. ' 

In contrast to common approaches in the literature, the pheromone update rule (18) does not 
consider the value of the objective function of a solution. In the distributed planning problem at hand, 
the agents do not reveal their individual costs to the mediator (cf. Section 2) who controls the search 
process and updates the pheromone values. Therefore, the pheromone update depends largely on the 
encoded solution jointly accepted by the agents. 



(18) 



4 Evaluation 

4.1 Setup of computational study 

4.1.1 Goals 

The metaheuristic CACM is evalauted by means of a computational benchmark study. Three goals are 
pursued. First, the impact of the three voting rules from Sect. 3.6 is studied. Second, the effect of the 
new collaborative local search heuristic CLS from Sect. 3.5 is tested. Third, the performance of CACM 
is compared to results from the literature, that is, to the best known solutions for the non-distributed 
MLULSP and other heuristics for the distributed MLULSP. 

4.1.2 Instances 

All in all, 528 benchmark instances are used in this study^. Veral and LaForge (1985), Coleman and McKnew 
(1991), and Dellaert and Jeunet (2000) introduced 176 instances for the non-distributed MLULSP. In 
these instances, a single agent (|^| = 1) is responsible to produce all items and the instances are di- 
vided into three groups of small, medium, and large instances denoted as si, ml, and 11 (cf. Table 2). 
Based on these MLULSP instances, Homberger (2010) introduced 176 instances were the items are 
jointly produced by two agents (|^| = 2) and 176 instances were the items are jointly produced by five 
agents {\A\ = 5). Accordingly, these 352 DMLULSP instances are divided in six groups denoted as s2, 
s5, m2, m5, 12, and 15 (cf. Table 2). In instance group s5, each agent produces the same number of 
items, while in the remaining groups (s2, m2, m5, 12, and 15) this is not the case. 

4.1.3 Measure of solution quality 

In the following, the objective function value f{y) (10) of the DMLULSP is considered in the context of 
the non-distributed MLULSP which is an approach suggested by Dudek and Stadtler (2005, 2007). For 
a given instance, the computed solution y for the distributed MLULSP is compared to the best-known 
solution y^'' for the non-distributed MLULSP and the percentage gap Q{y) is calculated as 

^Instances are available at http://www.dmlulsp.com. 
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Table 2: Characteristics of used benchmark instances 



class 


group 


1^1 

1 1 


m 


n 


no. of instances 


s 




1 
1 


o 


1 9 


aV 






9 


K 
O 


1 9 








o 


K 
O 


1 9 

X z 


w 


m 


ml 


1 


40, 50 


12, 24 


40 




m2 


2 


40, 50 


12, 24 


40 




m5 


5 


40, 50 


12, 24 


40 


1 


11 


1 


500 


36, 52 


40 




12 


2 


500 


36, 52 


40 




15 


5 


500 


36, 52 


40 



Table 3: Parameter during evaluation of metaheuristic CACM 





class s 


class m 


class 1 




1 


1 


1 


^max 


100 


1,000 


1,500 


P 


5% 


5% 


5% 


a 


5% 


5% 


5% 


s 


50,000 


200,000 


400,000 



^(y) = ^^';wfjlf'^ 100- (19) 

The best-known values f"''^{y^^) have been gathered from Dellaert and Jeunet (2000), Pitakaso et al. 
(2007), Homberger (2008, 2010), and Xiao et al. (2011). For the instances in si, the optimal solutions 
are known. 



4.1.4 Implementation and hardware 

CACM was implemented in JAVA (JDK 1.6). The computational experiments were executed on a 
Linux personal computer with Intel Core 2 Duo processor 1.83 GHz. To compare the heuristics in 
Sect. 4.4, all experiments were executed on this computer. 



4.1.5 Parameters 

By means of a preliminary study, the parameters of CACM shown in Table 3 have been determined. 
For all instances, the pheromone evaporation rate p, the pheromone intensification factor a, and the 
minimum pheromone value r™"*"" per edge are constant. Then again, the number of generated solutions 
s and the maximum pheromone value r*""^ per edge grow with increasing instance size. The ballot 
size b is evaluated together with the voting rules in the next section. 



4.2 Effect of voting rules and ballot size 

The effects of the Borda and Rawls voting rules on 
of the benchmark instances from groups s2 and s5. 
of 1, 10, 20, and 100. The results are given in Table 



the performance of CACM are studied by means 
Both voting rules are applied using ballot sizes b 
4. 
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Table 4: Effect of Borda and Rawls voting 



Agents Ballot size b Q 







Borda 


Rawls 


\A\ = 2 


1 


2.21 


2.21 




10 


2.14 


2.13 




20 


2.40 


2.31 




100 


2.46 


2.55 


1^1=5 


1 


8.31 


8.31 




10 


8.78 


8.98 




20 


8.62 


9.21 




100 


9.22 


9.41 



Table 5: Effect of collaborative local search. 



Heuristic 




G 




si 


s2 


s5 


CLS 


3.93 


5.29 


9.53 


CACM"^' 


1.65 


4.96 


9.15 


CACM 


0.07 


2.21 


8.31 



The configuration of Borda voting and Rawls voting with a ballot size of 6 = 1 results in equivalent 
voting procedures and, therefore, leads to identical results (see Table 4). Especially in case of scenarios 
with five agents, the results for the Borda voting rule are slightly better than the results for the Rawls 
voting rule. However, in the five agents scenario, both voting rules lead to the best results for ballot 
sizes of only one. For the two agent scenario, ballot sizes of one lead to the second best results. All in 
all, both voting rules seem to result in higher total costs compared to the degenerated case (6 = 1). 
Therefore, the parameter ballot size b is set to one, that is, the function 'voting(i?, A)' in Alg. 1 is in 
fact deactivated. The only voting mechanism applied is complete approval voting which is used during 
collaborative local search. A positive side effect of doing without Rawls and Borda voting is that no 
local cost information have to be revealed (Rawls) and the opportimity of strategic bidding, which is 
intrinsic in ranked based voting, is reduced (Borda). 

4.3 Effect of collaborative local search 

In this test, the effect of CLS is studied. The CACM approach including collaborative local search 
(CLS, cf. Sect. 3.5) is compared with CACM excluding CLS. The latter is denoted as CACM^^'. As 
described previously, only complete approval voting is used. Furthermore, CLS is used as stand-alone 
heuristic which follows Alg. 1 without the 'constructEncodedSolution' procedure. The initial encoded 
solution e is set to the unit matrix. 

The three variants are compared by means of the instances from groups si, s2, and s5. The 
instances of si correspond to the non-distributed MLULSP for which optimal solutions are known. 
Consequently, Table 5 depicts the average gap Q with respect to the optimal MLULSP solution. 

CACM solves 80 of 96 instances to optimality. For the remaining 16 instances, the average optimal- 
ity gap is only 0.07 percent although CACM does not use heuristic information specific to the MLULSP 
to control the search process. Furthermore, CACM outperforms the stand-alone CLS and CACM^^ 
both in the non-distributed (si) and distributed (s2, s5) scenario. Obviously, the combination of CLS 
and CACM^^ increases the overall performance. Finally, Figure 3 shows exemplarily for an instance, 
that the heuristic CACM converges faster than the heuristic CACM*^^. 
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Figure 3: Convergency of CACM and CACM' 



Table 6: Comparison of different heuristics by means of average percentage gap Q from best-known 
non-distributed solutions. 







Class s 






Class m 






Class 1 






ES 


SA 


AGO 


CACM 


ES 


SA 


ACQ 


CACM 


ES 


SA 


ACO 


CACM 


1^1 = 2 
|A| = 5 


5.92 
18.66 


1.78 
3.18 


2.66 
7.63 


2.21 
8.31 


18.65 
55.38 


7.95 
16.55 


2.77 
8.18 


1.84 
7.74 


55.07 
95.93 


16.84 
46.38 


13.30 
23.79 


5.43 
5.17 



4.4 Comparison with heuristics from the Hterature 

CACM is compared to three approaches designed for the DMLUSLP denoted here as SA, ES, and ACO 
which are based on the metaheuristics simulated annealing, evolutionary strategy, and ant colony opti- 
mization, respectively. SA was introduced by (Homberger, 2010), ES is presented in Homberger (2011), 
and ACO is discussed in (Homberger and Gehring, 2010, 2011). Some of these previous computational 
experiments were repeated or broadened to take into account additional instances (see Section 4.1.4). 

Table 6 shows the average gap Q over 176 instances with two agents (classes s2,m2,12) as well as over 
176 instances with five agents (classes s5, m5, 15). The detailed results for medium and large instances 
are given in Tables 7, 8, 9, and 10 in the appendix. Due to long runtimes Homberger and Gehring 
(2010, 2011) compared their approach only by means of the medium sized instances of groups m2 and 
m5. 

The presented CACM approach significantly improves the best known solutions for the medium 
and large instances (see Table 6). For example, CACM reduces the average gap Q of the large instances 
with five agents (group 15) from 46 percent to only about 5 percent. In detail, CACM improves twenty- 
nine out of forty m2 instances, twenty out of forty m5 instances, thirty-six out of forty 12 instances, 
and forty out of forty 15 instances (see appendix. Tables 7, 8, 9, and 10). With respect to the classes 
of small instances s2 and s5, CACM is the second best procedure. 
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5 Conclusions 



A collaborative solution approach based on the metaheuristic ant colony optimization was presented. 
It is used to coordinate joint production planning of a coalition of self-interested decision makers 
with asymmetric information. The search process is executed by a mediator who is not aware of the 
agents local costs. Furthermore, the agents do not have to reveal private cost information during 
search. Nevertheless, to control the search process, the mediator tries to elicit preferences by means of 
voting and a new collaborative local search procedure. The results are used to update the pheromone 
information which guide the solution construction process. Solutions are constructed by means of a new 
graph structure based on encoded solutions which is less complicated and therefore easier to search. 
The new approach CACM is evaluated by means of 356 benchmark instances from the literature and 
outperforms all existing distributed planning approaches on the sets of medium and large instances. 

For future research, it appears interesting to study collaborative lot sizing problems with a non- 
disjunct assignment of items to agents. This would increase potential planning conflicts between the 
agents and opens up interesting challenges for developing collaborative planning metaheuristics. 
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Appendix: Detailed Computational Results 



Table 7: Results per instance of class m with two agents {\A\ = 2) 
1^1 — 1 distributed to \A\ = 2 agents 



r*l TD i Tl Qt n Tl PA 


UCo L JVllU W 11 


ES 


SA 


AGO 


CACM 


1 


194,571.00 


226,310.55 


199,961.45 


202,824.35 


198,913.65 


2 


165,110.00 


189,359.10 


172,215.95 


168,960.35 


171,589.85 


3 


201,226.00 


243,635.40 


207,379.75 


207,244.65 


201,436.75 


4 


187,790.00 


224,011.65 


191,731.45 


191,036.65 


189,687.65 


5 


161,304.00 


184,229.85 


165,815.30 


166,007.00 


164,463.90 


6 


342,916.00 


397,948.15 


351,960.65 


349,865.05 


351,508.00 


7 


292,908.00 


338,136.05 


305,473.45 


299,409.35 


297,715.40 


8 


354,919.00 


422,567.70 


364,522.50 


367,615.55 


362,346.40 


9 


325,212.00 


357,269.80 


333,153.40 


334,031.75 


332,133.80 


10 


385,939.00 


441,387.35 


407,565.35 


398,898.65 


393,280.40 


11 


179,762.00 


192,527.35 


189,529.80 


204,292.60 


179,839.00 


12 


155,938.00 


164,848.85 


182,688.95 


157,301.95 


156,130.95 


13 


183,219.00 


191,549.55 


189,108.80 


188,364.60 


183,242.25 


14 


136,462.00 


145,731.45 


140,651.90 


138,738.35 


141,477.60 


15 


186,597.00 


214,290.70 


191,057.20 


190,761.45 


188,255.45 


16 


340,686.00 


397,568.85 


351,210.65 


357,909.70 


349,256.35 


17 


378,845.00 


448,356.00 


393,817.25 


391,014.60 


382,513.70 


18 


346,313.00 


411,172.95 


372,168.80 


357,764.35 


356,937.55 


19 


411,997.00 


464,478.90 


442,328.20 


421,823.70 


417,107.60 


20 


390,194.00 


426,386.00 


412,580.70 


402,375.25 


407,367.20 


21 


148,004.00 


214,663.10 


170,440.75 


148,628.90 


148,083.70 


22 


197,695.00 


248,894.10 


241,675.30 


201,297.00 


200,501.70 


23 


160,693.00 


201,806.10 


178,578.30 


160,924.90 


160,924.90 


24 


184,358.00 


216,154.35 


220,872.40 


185,426.70 


185,157.20 


25 


161,457.00 


205,387.90 


177,611.55 


163,421.00 


161,967.00 


26 


344,970.00 


448,962.10 


404,290.70 


349,062.25 


349,649.20 


27 


352,634.00 


471,776.05 


403,475.40 


363,633.60 


369,683.25 


28 


356,323.00 


480,489.55 


435,463.50 


363,952.10 


361,585.10 


29 


411,338.00 


651,772.35 


499,775.40 


431,998.00 


429,759.90 


30 


401 732 00 


521 972 05 


442 778 65 


41 4 74S 1 


414 996 80 


31 


185 161 00 


210 866 20 


200 677 75 


1 89 066 75 


191 642 55 


32 


185,542.00 


207,226.15 


193,634.60 


191,893.95 


188,634.85 


33 


192,157.00 


211,073.20 


204,163.95 


201,108.45 


195,519.95 


34 


136,757.00 


152,061.20 


147,227.50 


138,928.40 


138,987.85 


35 


166,041.00 


181,383.50 


177,690.45 


170,896.25 


168,874.65 


36 


289,846.00 


332,234.45 


308,642.60 


292,530.95 


293,472.75 


37 


337,913.00 


384,122.10 


354,226.60 


348,464.80 


344,700.05 


38 


319,905.00 


377,400.80 


342,019.30 


327,445.40 


328,214.05 


39 


366,848.00 


424,664.05 


387,213.45 


377,137.50 


374,957.45 


40 


305,011.00 


380,260.00 


319,977.20 


310,249.80 


308,190.20 


median 


245,536.00 


290,564.28 


273,574.38 


249,887.80 


247,454.75 


mean 


263,157.33 


315,123.39 


284,383.92 


270,676.37 


268,517.66 


stand, dev. 


94,057.58 


125,416.02 


104,568.95 


97,388.37 


97,485.57 


median Q 




15.90 


5.67 


2.55 


1.69 


mean Q 




18.65 


7.95 


2.77 


1.84 


stand, dev. Q 




10.64 


5.96 


2.09 


1.28 


no. of best known" 










11 


30 



''Both AGO and CACM compute the best known solution for instance 23. 
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Table 8: Results per instance of class m with five agents {\A\ = 5) 





\A\ = 1 




distributed to 


= 5 agents 




class m instance 


best-known 


ES 


SA 


AGO 


CACM 


1 


194,571.00 


264,065.00 


205,787.15 


204,753.85 


207,630.80 


2 


165,110.00 


231,090.80 


173,471.95 


173,361.15 


174,708.10 


3 


201,226.00 


280,535.50 


226,796.85 


214,696.05 


214,744.65 


4 


187,790.00 


273,135.80 


196,780.85 


205,980.65 


203,166.05 


5 


161,304.00 


200,282.10 


173,262.45 


174,641.15 


168,236.50 


6 


342,916.00 


493,419.20 


368,168.85 


369,897.35 


370,721.90 


7 


292,908.00 


403,347.70 


320,741.40 


316,507.00 


313,263.90 


8 


354,919.00 


563,164.80 


376,777.75 


373,085.55 


373,117.70 


9 


325,212.00 


447,899.60 


345,253.25 


359,006.95 


359,566.55 


10 


385,939.00 


551,828.05 


441,193.25 


425,803.55 


429,034.55 


11 


179,762.00 


312,657.85 


215,994.40 


200,473.80 


197,660.40 


12 


155,938.00 


198,585.80 


173,838.45 


171,895.75 


168,395.75 


13 


183,219.00 


251,756.45 


227,887.45 


207,910.65 


202,287.10 


14 


136,462.00 


165,146.70 


142,493.75 


148,623.95 


153,552.15 


15 


186,597.00 


343,134.10 


193,885.00 


203,231.10 


198,769.05 


16 


340,686.00 


541,487.35 


377,353.90 


377,807.85 


377,726.10 


17 


378,845.00 


596,551.85 


458,401.05 


413,040.65 


414,022.55 


18 


346,313.00 


525,989.70 


405,675.05 


391,587.65 


390,264.35 


19 


411,997.00 


559,655.00 


456,057.80 


468,349.70 


468,899.25 


20 


390,194.00 


600,109.20 


480,153.95 


457,261.15 


449,742.00 


21 


148,004.00 


211,273.15 


186,605.50 


154,271.80 


153,865.00 


22 


197,695.00 


324,509.20 


239,446.35 


205,484.50 


204,578.80 


23 


160,693.00 


261,299.90 


178,689.70 


171,443.70 


168,981.65 


24 


184,358.00 


303,488.65 


226,545.15 


190,604.80 


188,964.45 


25 


161,457.00 


248,752.95 


210,579.05 


169,652.40 


168,805.40 




■iAA Q7n nn 


77c; Kofi qc 


4'H4 79"? fin 


'H7fi ■^9'i Qn 


OLIO, 1 J-^.OL/ 


27 


352,634.00 


768,345.20 


596,826.80 


385,677.45 


383,341.15 


28 


356,323.00 


803,754.00 


464,721.50 


380,075.70 


383,552.95 


29 


411,338.00 


845,876.45 


718,175.10 


459,726.25 


459,557.10 


30 


401,732.00 


833,605.10 


458,301.15 


432,181.35 


429,949.70 


31 


185,161.00 


249,657.65 


207,223.10 


201,043.10 


198,819.55 


32 


185,542.00 


262,525.70 


205,138.90 


195,533.20 


197,363.30 


33 


192,157.00 


Otto 100 OA 

258,188.80 


Oil A/^ A OA 


OAC OOC; CA 

205,325.50 




■^4 


1 'Hfi 7^17 nn 

±OU, 1 O 1 .uu 


91 n QSfi 1 








35 


166,041.00 


238,389.35 


178,082.85 


179,720.25 


181,454.00 


36 


289,846.00 


419,973.85 


339,317.90 


308,004.50 


309,195.95 


37 


337,913.00 


464,681.45 


372,082.80 


359,501.40 


359,770.00 


38 


319,905.00 


443,159.95 


343,042.05 


342,613.85 


342,809.45 


39 


366,848.00 


534,803.25 


418,399.85 


397,048.95 


393,949.00 


40 


305,011.00 


498,094.70 


328,389.55 


321,507.00 


317,653.55 


median 


245,536.00 


373,240.90 


280,093.88 


261,350.27 


261,970.30 


mean 


263,157.33 


419,019.72 


310,837.55 


285,878.95 


284,825.75 


stand, dev. 


94,057.58 


194,238.04 


133,327.52 


105,938.42 


105,897.62 


median G 




45.17 


11.70 


8.14 


7.09 


mean Q 




55.38 


16.55 


8.18 


7.74 


stand, dev. Q 




26.49 


14.65 


3.02 


2.92 


no. of best known 







8 


12 


20 



18 



Table 9: Results per instance of class I with two agents (|^| = 2) 





\A\ = 1 




distributed to \A\ =2 agents 




class / instance 


best-known 


ES 


SA 


AGO 


CACM 


1 


591,585.00 


695,337.65 


608,572.30 


642,747.80 


597,566.10 


2 


816,043.00 


923,305.30 


885,146.00 


878,505.40 


982,418.75 


3 


908,616.00 


1,113,228.95 


975,174.80 


1,000,429.40 


996,370.50 


4 


929,929.00 


1,037,113.85 


1,035,097.40 


1,030,539.85 


1,011,923.30 


5 


1,145,749.00 


1,371,152.00 


1,255,854.25 


1,222,565.20 


1,266,812.85 


6 


7,639,920.00 


12,781,176.60 


9,782,048.90 


8,510,346.95 


8,198,323.90 


7 


3,921,407.00 


5,785,793.10 


5,275,859.90 


4,327,733.60 


4,054,305.25 


8 


2,694,924.00 


3,586,973.70 


3,235,243.15 


2,870,682.50 


2,824,140.80 


9 


1,880,021.00 


2,453,164.65 


2,061,679.05 


2,128,470.80 


1,899,801.85 


10 


1,502,194.00 


1,901,132.75 


1,624,990.25 


1,691,979.45 


1,569,679.65 


11 


59,121,845.00 


120,529,283.50 


73,444,728.85 


64,597,231.15 


59,407,117.90 


12 


13,422,827.00 


22,643,066.50 


15,480,871.65 


14,510,728.55 


13,573,477.20 


13 


4,718,146.00 


7,363,067.15 


5,867,821.10 


5,246,854.30 


5,008,244.80 


14 


2,908,634.00 


3,819,216.00 


3,195,148.95 


3,202,206.60 


3,068,236.60 


15 


1,737,525.00 


2,198,085.20 


1,853,347.30 


1,962,320.40 


1,838,495.35 


16 


468,463,630.00 


959,595,328.50 


555,959,053.05 


520,424,836.30 


473,526,552.45 


17 


18,677,678.00 


41,742,064.65 


21,357,104.50 


20,234,057.05 


18,860,888.70 


18 


7,308,193.00 


12,586,168.40 


8,640,468.25 


7,917,765.70 


7,526,011.90 


19 


3,519,932.00 


4,965,055.40 


4,028,118.60 


3,947,226.30 


3,713,975.95 


20 


2,278,214.00 


2,803,612.00 


2,446,637.20 


2,549,066.10 


2,401,407.05 


21 


1,187,090.00 


1,407,774.55 


1,323,941.05 


1,303,687.95 


1,243,593.70 


22 


1,341,584.00 


1,559,006.70 


1,474,146.65 


1,493,191.95 


1,536,134.45 


23 


1,400,480.00 


1,687,215.15 


1,518,462.70 


1,548,463.70 


1,502,495.80 


24 


1,382,150.00 


1,639,929.55 


1,533,006.85 


1,563,371.50 


1,499,628.70 


25 


1,657,248.00 


1,986,023.70 


1,939,231.90 


1,859,094.20 


1,814,371.90 




1 9 (S71 sns nn 


9c: 1 0K c^^q on 


15,998,629.30 


14,568,749.65 


1 H 270 'lO'i 1 


27 


7,159,416.00 


12,131,464.00 


8,501,797.65 


8,461,161.55 


7,496,052.05 


28 


4,148,783.00 


5,792,328.05 


A 00/1 A10 CK 

4,824,Uiz.00 


4, / /U, 983.10 


4,423,193.05 


29 


2,889,151.00 


3,776,572.00 


3,117,732.50 


3,339,847.50 


3,020,055.65 


30 


2,183,815.00 


2,732,881.40 


2,378,237.60 


2,588,377.55 


2,315,291.70 


31 


101,497,679.00 


278,824,944.75 


128,549,165.95 


127,935,080.15 


102,255,225.00 


32 


18,028,225.00 


37,609,906.05 


20,189,757.75 


21,578,019.85 


18,287,579.60 


33 


6,780,986.00 


11,353,745.60 


1 n rino TAG on 

10,902,146.20 


Apfi 010 on 

8,401,813.80 


7,318,975.85 


34 


4,055,536.00 


5,341,405.15 


5,029,223.90 


4,927,903.20 


4,283,010.35 


35 


2,559,885.00 


3,307,932.90 


2,810,070.40 


3,050,680.85 


2,771,436.85 


36 


755,506,278.00 


2,036,463,721.10 


955,384,033.75 


863,134,800.00 


768,264,766.55 


37 


33,309,777.00 


86,653,413.25 


45,203,262.65 


37,584,305.00 


33,724,916.75 


38 


10,464,662.00 


19,257,108.00 


12,772,705.35 


11,915,138.85 


10,734,284.00 


39 


5,116,338.00 


7,422,470.10 


6,072,238.50 


6,041,357.80 


5,441,242.45 


40 


3,391,440.00 


4,574,457.65 


3,666,782.10 


4,133,474.95 


3,515,181.75 


median 


3,455,686.00 


4,769,756.53 


3,847,450.35 


4,040,350.63 


3,614,578.85 


mean 


39,522,983.58 


93,963,529.39 


48,805,038.77 


44,977,394.91 


40,176,092.25 


stand, dev. 


136,399,204.84 


347,139,859.03 


169,755,185.52 


154,825,527.76 


138,439,704.88 


median Q 




33.99 


14.39 


12.16 


5.10 


mean Q 




55.07 


16.84 


13.30 


5.43 


stand, dev. G 




43.75 


10.62 


4.72 


3.96 


no. of best known 







2 


2 


36 



19 



Table 10: Results per instance of class I with five agents {\A\ = 5) 





\A\ = 1 




distributed to 


yl| = 5 agents 




class I instance 


best-known 


ES 


SA 


AGO 


CACM 


1 


591,585.00 


804,893.65 


710,850.55 


647,247.70 


635,739.60 


2 


816,043.00 


1,147,738.65 


988,143.20 


930,055.95 


821,528.70 


3 


908,616.00 


1,355,186.90 


1,188,264.05 


1,084,152.75 


948,842.40 


4 


929,929.00 


1,247,533.15 


1,184,439.70 


1,128,669.65 


967,533.50 


5 


1,145,749.00 


1,413,376.55 


1,259,973.25 


1,399,061.45 


1,205,951.95 


6 


7,639,920.00 


25,878,434.65 


19,037,568.40 


8,694,605.45 


8,355,460.70 


7 


3,921,407.00 


9,132,315.80 


8,707,302.60 


4,653,392.50 


4,237,945.50 


8 


2,694,924.00 


4,358,121.20 


3,438,619.75 


3,346,616.65 


2,947,906.85 


9 


1,880,021.00 


2,710,619.00 


2,137,343.75 


2,336,763.45 


2,028,723.05 


10 


1,502,194.00 


2,085,083.65 


1,677,771.85 


1,867,070.70 


1,583,783.40 


11 


59,121,845.00 


137,258,807.55 


98,380,736.30 


74,420,057.30 


59,486,366.90 


12 


13,422,827.00 


31,575,365.70 


18,208,444.20 


17,050,429.70 


13,619,454.45 


13 


4,718,146.00 


8,582,950.70 


6,387,892.00 


5,742,169.60 


5,117,738.80 


14 


2,908,634.00 


4,513,204.40 


3,650,417.90 


3,651,103.30 


3,172,980.55 


15 


1,737,525.00 


2,549,510.10 


1,951,706.95 


2,266,967.25 


1,828,004.20 


16 


468,463,630.00 


1,173,231,850.40 


724,511,376.10 


556,757,116.65 


473,495,284.60 


17 


18,677,678.00 


48,122,890.00 


31,079,449.25 


20,942,784.70 


18,866,385.50 


18 


7,308,193.00 


15,131,195.70 


10,952,393.35 


8,489,683.80 


7,641,505.35 


19 


3,519,932.00 


5,773,626.55 


4,671,451.50 


4,434,413.15 


3,811,046.35 


20 


2,278,214.00 


3,433,568.70 


2,855,260.20 


3,021,621.50 


2,466,302.90 


21 


1,187,090.00 


1,789,713.70 


1,618,514.00 


1,393,015.70 


1,232,991.95 


22 


1,341,584.00 


1,970,678.35 


1,660,601.70 


1,637,945.90 


1,400,411.50 


23 


1,400,480.00 


1,902,855.80 


1,694,771.35 


1,666,318.55 


1,416,472.45 


24 


1,382,150.00 


1,889,928.60 


1,583,299.20 


1,780,922.50 


1,406,540.30 


25 


1,657,248.00 


2,103,313.15 


1,987,959.25 


2,070,126.30 


1,739,024.90 




1 9 fi7i sns nn 




oc; QRq yon 0'^ 


1 fi 071 7Q'\ 8^ 




27 


7,159,416.00 


15,489,176.70 


12,806,323.95 


9,074,157.40 


7,780,993.10 


28 


4,148,783.00 


7,384,138.30 


5,620,887.95 


5,449,223.35 


4,657,224.70 


29 


2,889,151.00 


4,480,227.95 


3,571,646.00 


3,882,516.70 


3,011,682.40 


30 


2,183,815.00 


3,160,428.90 


2,567,987.70 


2,859,877.00 


2,264,261.75 


31 


101,497,679.00 


417,268,369.80 


302,711,993.50 


128,419,919.00 


102,344,312.65 


32 


18,028,225.00 


58,285,339.85 


35,763,440.10 


23,474,329.35 


18,573,788.45 


33 


6,780,986.00 


14,718,304.85 


11,170,886.30 


O Af\A Otto AA 

8,494,852.00 


7, 396,134. 75 




4:,UO?J,(JOU.UU 


1 ,U 1 fJ,UOU.OfJ 




0,OifJ,^O^.ifJ 


4 ^44 sn^ 

1 , •ill .('V'r> ••ill 


35 


2,559,885.00 


4,043,062.15 


2,968,750.20 


3,475,333.80 


2,680,782.55 


36 


755,506,278.00 


2,594,443,728.80 


1,259,991,223.40 


863,673,114.45 


768,319,555.25 


37 


33,309,777.00 


101,154,724.65 


51,080,825.70 


37,632,659.60 


33,725,180.95 


38 


10,464,662.00 


24,143,971.55 


15,436,884.20 


12,386,857.75 


10,804,778.05 


39 


5,116,338.00 


9,116,234.95 


8,180,743.30 


6,700,266.75 


5,678,066.85 


40 


3,391,440.00 


5,232,692.40 


3,935,280.35 


4,556,441.20 


3,557,098.90 


median 


3,455,686.00 


5,503,159.48 


4,303,365.93 


4,495,427.18 


3,684,072.63 


mean 


39,522,983.58 


119,757,107.01 


67,454,377.75 


46,573,473.48 


40,221,927.43 


stand, dev. 


136,399,204.84 


440,099,405.84 


225,778,836.85 


157,709,294.72 


138,436,787.71 


median Q 




62.87 


31.75 


24.60 


4.80 


mean Q 




95.93 


46.38 


23.79 


5.17 


stand, dev. Q 




69.34 


39.91 


6.76 


3.15 


no. of best known 













40 



20 



